Screening of magnetic moments in PuAm alloy : LDA+DMFT study 



J. H. Shim 1 , K. Haulc 1 , S. Savrasov 2 , and G. Kotliar 1 

1 Department of Physics, Rutgers University, Piscataway, NJ 08854, USA and 
2 Department of Physics, University of California, Davis, CA 95616 
(Dated: June 2, 2008) 

The puzzling absence of Pu magnetic moments in a PuAm environment is explored using the 
self-consistent Dynamical Mean Field Theory (DMFT) calculations in combination with the Local 
Density Approximation. We argue that 5-Pu -Am alloys provide an ideal test bed for investigating 
the screening of moments from the single impurity limit to the dense limit. Several important effects 
can be studied: volume expansion, shift of the bare Pu on-site / energy level, and the reduction of 
the hybridization cloud resulting from the collective character of the Kondo effect in the Anderson 
lattice. These effects compensate each other and result in a coherence scale, which is independent 
of alloy composition, and is around 800 K. We emphasize the role of the DMFT self-consistency 
condition, and multiplet splittings in Pu and Am atoms, in order to capture the correct value of the 
coherence scale in the alloy. 



PACS numbers: 

The electronic structure of plutonium (Pu) is one of the 
outstanding issues in condensed matter theory. The vol- 
ume of <5-Pu is between that of the early actinides, where 
the / electrons are itinerant, and the late actinides, where 
the / electrons are localized. Hence the view that the / 
electrons are very close to a localization-delocalization 
transition was put forward early on by Johansson^. 

Theoretical density functional theory calculations in 
many implementations, all predict <5-Pu to be magnetic^. 
Nevertheless, various experiments that searched for the 
magnetic moment in <5-Pu, such as specific heat in field, 
neutron scattering and /zSR measurements^, found no 
evidence for the ordered, or fluctuating moment in the 
energy and temperature windows probed by these exper- 
iments. High energy spectroscopics, such as electron- 
energy-loss spectroscopy and X-ray absorption spec- 
troscopy experiments of the 5d-5/ core-valence transi- 
tion, clearly indicate that / electrons in Pu are in a 
5/ 5 configuration^, which carries a magnetic moment. 
Reconciling these low energy and high energy experi- 
ments, has lead to numerous conflicting theories. Within 
the most recent Dynamical Mean Field Theory (DMFT) 
studies of <5-Pu, allowing for magnetic and non magnetic 
solutions, Pu does not order magnetically, and instead 
forms a non-magnetic Fermi liquid with a coherence tem- 
perature T* of the order of 800 K^. The Pu atoms are 
in a mixed valent state, with a dominant 5/ 5 configu- 
ration, but with sizeable admixture of 5/ 6 . On the ba- 
sis of this picture it was expected, and shown explicitly 
via a landmark continuous-time quantum Monte Carlo 
calculation^, that expanding the lattice would reduce the 
coherence temperature T* , making the moment of the / 5 
configuration observable. Other theories, whereby the / 
moments are quenched directly by spin pairing, which is 
mediated by the spin orbit coupling, also account for the 
absence of moments in Pu, and suggest that this material 
is at a quantum critical point£. 

PuAm alloy is a unique system to study the effect of 
lattice expansion of <5-Pu since the volume of Amcricium 
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FIG. 1: The volume of Pu and Am in each phases. According 
to Ref.[^, the volume of Puo.76 Amo.24 is locates slightly above 
(~ 2%) the value estimated from the linear interpolation be- 
tween 5-Pu and Am-I. 



at ambient pressure (Am-I) is 20% bigger than volume 
of <5— Pu (see Fig. [I]). Moreover, there is no structural 
phase transition between pure 5— Pu and Puo.25Amo.75, 
and the structure remains /coi^. Surprisingly, an in- 
tensive program studying the Pu-Am mixtures, which 
has the effect of expanding the distance between the Pu 
atoms, has not shown an enhanced magnetic susceptibil- 
ity, nor a narrowing of the peak at the Fermi level, as 
measured by photoemission experiment&ii. This seems 
naively at odds with both the quantum critical picture, 
and the DMFT picture, and supports the suggestion, that 
the configuration of Pu in the solid is close to an inert 
5/ 6 configuratio n 1 2 1 13 , being essentially equivalent to Am 
configuration. Unfortunately, this picture can not ex- 
plain the results of high energy spectroscopies^, nor the 
difference in specific heat between ev-Pu and <5-Pu&^, 
nor the spectral properties^ 5 .. Hence the puzzle of the 
absence of magnetic moments in Pu-Am alloys remains. 

In this paper, we provide a theoretical explanation of 
this puzzle using LDA+DMFT methodic. We find 
that indeed the moment remains well screened in agree- 
ment with experiment, and we identify the mechanism 
that compensates the effects of the lattice expansion. 
While the Pu-Am system cannot realize the unscreening 
of the local moment, and hence the approach to quantum 
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criticality, it provides an ideal system in which to follow 
the transition from collective screening of moments, to 
single impurity spin screening of moments. 

Within the LDA+DMFT method^, the itinerant spd 
electrons are treated by the LDA method^, while for 
the electrons in the correlated / orbitals, all the local 
Feynman diagrams are summed up by solving an aux- 
iliary quantum impurity problem, subject to a DMFT 
self-consistency condition 
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Here A(w) is the hybridization matrix of the auxiliary 
impurity problem, Ei mp is the matrix of impurity lev- 
els, and Gk is the one electron Green's function G*k = 
l/(Ok(w+//)-if k -£). -ffk and O k arc the LDA Hamil- 
tonian and the overlap matrix in the projective orthogo- 
nalized LMTO basis sefc^. The quantum impurity solver 
is used to determine the local self-energy correction to the 
/-orbital Hamiltonian, which is a functional of the impu- 
rity levels, and the hybridization strength Y\Ei mp , A]. To 
solve the impurity problem, we used the vertex corrected 
one-crossing approximation method 2 ^, and the results are 
further cross-checked by the continuous time quantum 
Monte Carlo method 2 ^. The Slater integrals F k (k = 2, 
4, 6) are computed by Cowan's atomic Hartree-Fock pro- 
gram including rclativistic corrections 2 ^ and rescaled to 
80% of their atomic value, accounting for the screening 
in the solid. We use Coulomb interaction U = 4.5 eV, 
consistent with previous studies of Pu and A m 23 ' 24 . 

The Pui-jjAmj alloys are approximated by the charge 
ordered compounds with x = 0, 1/4, 1/2, and 3/4. Based 
on the original fee phase, we use a four times larger unit 
cell with the simple cubic structure. To determine the 
lattice constant of the alloy, we used the linear interpola- 
tion between the two end compounds, 5-Yw and /3-Am^. 
We use 8x8x8 momentum points in the first Brillouin 
zone in combination with the tetrahedron method. 

Figure[2] (a-c) shows the calculated spectral function of 
Pui-jAnij;. The theoretical results are compared to the 
experimental spectra of <5-Pu alloyed with Arn^. Naively, 
one would expect that alloys of Pu and Am are magnetic, 
because the volume of the alloy is bigger than the volume 
of the elemental plutonium. Increasing the volume de- 
creases the hybridization between the 5/-orbital and the 
conduction spd-orbitals, and consequentially the Kondo 
coupling J should be reduced. Since the strength of the 
Kondo screening depends exponentially on the Kondo 
coupling J while RKKY scale is only a quadratic function 
of J, one is expecting the screening to become inefficient 
even for a small change of J. However, experiments ruled 
out the possibility of magnetic moment formation. The 
specific heat of Pui-^Am^. does not show any significant 
changes 1 * with alloying and is simply proportional to the 
Pu content. This is because the inert f 6 configuration of 
Am has very small specific heat, and only Pu atoms con- 
tribute to c v . The photocmission spectral of Pui-^Am^ 
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FIG. 2: (a-e) The spectral function of Pui-iAni^ in the fee 
unit cell. Total (black), Pu 5/ (red), Am 5/ (blue) spectrum 
are obtained by the LDA+DMFT method. The green lines 
show experimental photoemission spectrum at x = 0.0 (a), 
0.20 (b), 0.36 (c), and 1.0 (e) (extracted from Ref. 0). (f) 
The spectral function of Pu atom for various Am dopings. 



has been reproduced by the simple summation of the 
spectrum of pure <5-Pu and Am^. 

Our LDA+DMFT results in Fig. H show a very favor- 
able agreement with experiment. The spectral function 
of Pu atom, displayed in Fig. [3(f), is almost unchanged 
in PuAm alloys. The coherence scale, that determine the 
temperature below which the Pu moment is screened, re- 
mains roughly at 800 K in the alloys. The Am spectra 
is even more doping independent and shows no appreci- 
ate change when alloyed with Pu. This is not surprising 
since a very large pressure is needed to delocalize Am 
5/ electron a 24 ' 25 , and only the Am IV phase is known to 
be itinerant (see Fig.l for location of Am IV in volume 
diagram), corresponding to volume of Am IV in Fig. [TJ 
Taking into account that the spectra of both Pu and Am 
are almost unchanged in alloys, the total DOS can be 
well approximated by a linear combination of the DOS 
of the two end compounds, as conjectured in Ref. 9. 

In order to gain further insight into the puzzle of 
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the doping independent coherence scale in Pui^Aru^ 
compound, we will examine the important factors that 
govern the screening of the magnetic moments in our 
LDA+DMFT calculation. The temperature dependence 
of the spectra and hence the coherence scale is entirely 
determined by the DMFT self-energy correction E, which 
is a functional of the hybridization function A, impurity 
levels Ei mp and the screened Coulomb interaction U, i.e., 
£[A, Ei mp , U]. U is here kept constant. The hybridiza- 
tion function A and impurity levels E imp , determined by 
Eq. [TJ arc both quite sensitive to the volume change as 
well as chemical substitutions, as we will show below. 

The coherence scale can be expressed analytically when 
the following two assumptions are made: a) The atom 
has SU(N) symmetry and hence the spin-orbit, crystal 
field splittings and Hund's coupling are absent, b) the 
hybridization function A is frequency and temperature 
independent. In this case, the Kondo coherence scale 
can be estimated by T K ~ e -1 ' J , where J is 

n \E„+i - E n E n - E n _i) ^ ^ 

and E n are the atomic energies corresponding to occu- 
pancy n, N is degeneracy of the SU(N) model. In the ex- 
treme Kondo limit, the difference in the atomic energies 
arc E n +i — E n = U/2 and E n — E n -x = U/2. The above 
assumptions are not realistic in Pu because of the strong 
multiplet effects and the strong frequency and temper- 
ature dependence of the hybridization function A(u>,T). 
Nevertheless, the formula Eq. Ogives the hint that the co- 
herence scale is sensitive to the ratio (|A"|)/A.E/, where 
(| A" |) is average of the hybridization function in some 
low energy region, and AEf is the difference between 
the ground state of the atomic / 6 configuration and the 
ground state of the atomic / 5 configuration. Namely, 
Pu 5/ electrons are in mixed-valence state with an aver- 
age occupancy close to n/ w 5.3^ therefore the dominant 
fluctuations arc 
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FIG. 3: (a) The number of 5/ electrons, rif, as a function 
of volume (red dots), doping x at constant volume (green 
dots), and combination of both (blue dots). Note that a — 
8.1, 8.7 and 9.3 a.u. correspond to the volume of a-Pu, 5-Pu 
and /3-Am, respectively, (b) Volume dependence of hybridiza- 
tion function for Pu in fee phase, (c) The energy difference 
between the two most relevant atomic states, defined in the 
text. 



e + 5/ 5 ^5f 6 , (3) 

where e is a conduction spd electron. 

In turn, the 5/ occupancy strongly cooperates with the 
coherence scale, and the scale is smallest in the Kondo 
regime where rif is integer. In Fig. [3] (a) we show the 5/ 
occupancy for three different cases of chemical substitu- 
tion and volume change. If the volume of the elemental 
Pu is changed without alloying, n/ is slightly decreasing 
with increasing volume (see red curve in Fig. [3] (a)). This 
is to be expected since in the extreme limit of large vol- 
ume the system will be in the extreme Kondo limit with 
rif = 5. The magnitude of the rif variation in the range 
of volumes, corresponding to a-Pu and Am is however 
surprisingly small, only of the order of 0.06. If the vol- 
ume is kept constant, but the the chemical composition 
is varied, the occupancy rif is increasing with Am doping 
(green dots in Fig. [3J (a) correspond to x = 0, 1/4, 1/2, 
and 3/4 from left to right). The Pu atom is thus more 



mixed valent in the limit of small Pu concentration, i.e., 
single impurity limit. This might seem counter-intuitive, 
but the increase in the coherence scale can be understood 
from the following consideration. The quasiparticle peak 
detected in 5/ partial DOS creates a depletion in the 
conduction spd bands. This is traditionally called Kondo 
hole, which is always pinned at the Fermi level. The hy- 
bridization is thus decreased due appearance of quasipar- 
ticle peak, an effect termed by Noziere as "exhaustion" . 
This Kondo hole does not form on Am atom, since there 
is no coherence peak on Am atom. Thus, if a single Pu 
atom in emerged in a lattice of Am atoms, the moment 
can be more efficiently screened by spd-electrons, resid- 
ing on neighboring Am sites. Finally the blue dots in 
Fig. [3J (a) show that the two effects largely cancel, when 
both the substitution and volume change are taken into 
account, and the 5/ occupancy of Pui-a-Am^ is almost 
constant. 
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The frequency dependent hybridization function 
A"(cj)|/7t for elemental Pu at different volumes (no al- 
loying), is shown in Fig. [3] (b). The Kondo hole, i.e., 
depletion of hybridization at zero frequency, is clearly vis- 
ible. Moreover, the hybridization changes substantially 
when volume is varied between the volume of a-Pu and 
Am volume. The change is between 30% at low energy 
and up to 50% at intermediate energies. This is a sub- 
stantial change, which indeed unscrcens magnetic mo- 
ments and changes Pauli susceptibility into Curie- Weiss 
susceptibility, as shown in Ref. 0. 

The Eq. [5] reveals that the coherence scale depends 
sensitively on the ratio \A"\/ AEf rather than on the hy- 
bridization itself. In Fig. [3] (c) we plot the smallest energy 
difference AEf = Eq — £5, which is the dominant contri- 
bution to the exchange coupling in Eq.[5] Here Eq and £5 
are the ground state energies of the atomic / 6 and atomic 
/ 5 configuration, respectively. In Fig. [3] (c) we separate 
the effect of volume increase from the chemical substitu- 
tion. It is clear from Fig. [3] (c) that the energy differ- 
ence is almost three times smaller than U/2 ~ 2.25eV, 
the value that corresponds to the extreme Kondo limit. 
The 5/ 6 atomic configuration has thus substantial over- 
lap with the ground state wave function of the solid, as 
shown by valence histogram in Ref. @. It is also apparent 
from Fig. [3] (c) that the volume increase has only a small 
effect on this energy difference (red curve). When alloy- 
ing with Am, the change of AEf is larger (blue dots) and 
is of the order of 25%. Since AEf is decreasing, it thus 
partially compensates the decrease of hybridization with 
pressure, shown in Fig. 0(b). However, the magnitude of 
AEf variation is too small to compensate the dramatic 
change of hybridization. 
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FIG. 4: DMFT hybridization function for Pui-mAm^ at fixed 
volume (a) computed within the self-consistent LDA+DMFT 
method, (b) within LDA by setting E = 0. 



Finally, Fig. [4^ shows the hybridization function for 
each composition of Pui_ x Am x alloy for the case of a 
fixed volume. Because the dominant contribution to hy- 
bridization function is not the /-/ hybridization but the 
f-spd hybridization^, the overall magnitude of the hy- 
bridization is not changed. However, it is apparent from 
Fig. [4] (a) that the chemical substitution of Pu by Am 
eliminates the hybridization dip, and therefore the Pu 
moment is even more efficiently screened by the spd elec- 
trons from neighboring Am atoms than neighboring Pu 
atoms. This elimination of Kondo hole by the chemi- 
cal substitution by Am thus compensates the decrease of 
hybridization due to volume change. The net change of 
AEf and the change of hybridization (|A"|) due to both 
volume decrease and alloying compensate and lead to a 
constant coherence scale across the phase diagram. 

Fig. 5] (b) shows that the increase of hybridization due 
to alloying is present even on the LDA level, i.e., comput- 
ing A by Eq.Q]and neglecting the self-energy corrections. 
However, the magnitude of the hybridization A and the 
magnitude of the change of A is substantially larger in 
LDA than in self-consistent LDA+DMFT method. The 
"need for renormalization" of the LDA hybridization was 
pointed out very early o n 26 ' 27 in the context of the de- 
scription of cerium photocmission, when compared to a 
single impurity calculations. 

To conclude, LDA+DMFT accounts for all the salient 
features of the experiments, and give new insights into 
the mechanism for spin compensation in the late ac- 
tinides. It isolates two essential elements for the screen- 
ing of the magnetic moment: hybridization and position 
of the / level. While hybridization is decreased by vol- 
ume expansion, as naively expected, the substitution of 
Am eliminates the dip in the hybridization function, and 
at the same time shifts the position of / level. The three 
effects compensate each other, and the coherence scale 
remains around 800 K. Hence the Pu-Am mixtures ac- 
centuates the mixed valent character of Pu. 

Our work has important experimental consequences. 
While thermodynamic quantities such as susceptibility 
and specific heat should not exhibit significant modifica- 
tions upon alloying, the same is not true for transport 
quantities. Optical conductivity is an excellent probe of 
the hybridization gap, and we predict that upon alloying, 
this hybridization gap will diminish while the dc resistiv- 
ity will rapidly increases. 
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